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Abstract We describe a calculation of the spectrum of strange and nonstrange hadrons that simul- 
taneously correlates the dressed-quark-core masses of meson and baryon ground- and excited-states 
within a single framework. The foundation for this analysis is a symmetry-preserving Dyson-Schwinger 
equation treatment of a vector x vector contact interaction. Our results exemplify and highlight the 
deep impact of dynamical chiral symmetry breaking on the hadron spectrum: an accurate description 
of the meson spectrum entails a similarly successful prediction of the spectrum of baryons, including 
those with strangeness. The analysis also provides numerous insights into baryon structure. For exam- 
ple, that baryon structure is largely flavour-blind, the first radial excitation of ground-state baryons 
is constituted almost entirely from axial- vector diquark correlations, and DCSB is the foundation for 
the ordering of low-lying baryon levels; viz., (1/2)+, (1/2)+, (l/2)~. 

Keywords Bethe-Salpeter equation • Confinement • Dynamical chiral symmetry breaking • Dyson- 
Schwinger equations • Faddeev equation • Hadron spectrum • Strange mesons and baryons 



1 Introduction 

Quantum chromodynamics (QCD) is a relativistic quantum field theory that is generally believed 
to describe the strongly interacting part of the Standard Model. If so, it is Nature's only known 
example of an essentially nonperturbative fundamental theory. This is the difficulty: in attempting to 
understand QCD one must immediately confront a unique nonperturbative problem. Never before have 
we been confronted by a theory whose elementary excitations are not those degrees-of-freedom readily 
accessible via experiment; i.e., whose elementary excitations are confined. Moreover, it appears that 
QCD generates forces which are so strong that less-than 2% of a nucleon's mass can be attributed to the 
current-quark masses that appear in the QCD Lagrangian; viz., forces that generate mass from (almost) 
nothing, a phenomenon known as dynamical chiral symmetry breaking (DCSB). It follows that the 
Higgs mechanism is largely irrelevant to the bulk of normal matter in the Universe. Instead the most 
important mass generating mechanism is the strong interaction effect of DCSB. Neither confinement 
nor DCSB is apparent in QCD's Lagrangian and yet they play the dominant role in determining the 
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observable characteristics of real- world QCD. The physics of hadrons is ruled by emergent phenomena, 
such as these, which can only be elucidated through the employment of nonperturbative methods in 
quantum field theoryQ This is both the greatest novelty and the biggest challenge within the Standard 
Model. 

One method by which to validate QCD is computation of its hadron spectrum and subsequent 
comparison with modern experiment. Indeed, this is an integral part of the international effort in 
nuclear physics. The N* programme 0; Q and the search for hybrid and exotic mesons [3; O together 
address the questions: which hadron states and resonances are produced by QCD, and how are they 
constituted? Herein, motivated by this intense effort in hadron spectroscopy, we extend Rcf. Q and 
treat ground- and excited-state hadrons with s-quark content. Furthermore, as Ref. Q was a precursor 
to a wide-ranging study of nuclcon elastic and transition form factors 0; the study we describe 
herein is also a necessary step toward a comprehensive analysis of form factors that involve hadrons 
with strangeness. 

We br ing to these tasks a continuum perspective based on QCD's Dyson-Schwinger equations 
(DSEs) [11; IgUlOj and within this framework we use a symmetry-preserving treatment of a vector x vector 
contact interaction because that has proven to be a reliable tool in spectrum calculations. It is appro- 
priate to remark that this interaction produces form factors which are too hard 0;[l;[iil[i3 but, when 
interpreted carefully, they, too, can be used to draw valuable insights. 

To explain our choice of interaction we note by contrast that DSE kernels with a closer connection 
to perturbativc QCD; namely, which preserve QCD's one-loop rcnormalisation grou p b ehaviour, have 
long been employed in studies of the spectrum and interactions of mesons [l^ [l3; Il5| . Such kernels 
are developed in the rainbow-ladder approximation, which is the leading-order in a systematic and 
symmetry-preserving truncation scheme [TgI : [iTj ; and their model input is expressed via a statement 
about the nature of the gap equation's kernel at infrared momenta. With a single parameter that 
expresses a confinement length-scale or strength [l^ , they have successfully described and predicted 
numerous properties of vector [l^; [^; 21; 22; 23] and pseudoscalar mesons 0; [H; |2^; [U; [25l; [2^; [13] 
with masses less than 1 GeV, and ground-state baryons [1^; [1^ [S^; [Ml . Such kernels are also reliable 
for ground-state heavy-heavy mesons [3^ . 

On the other hand, whilst mo del- independent results for properties of pseudoscalar meson excited 
states have been established (SSl : |34| , the rainbow-ladder truncation is quantitatively inaccurate for the 
spectrum of light-quark mesons with masses greater than 1 GeV, for reasons which are understood [2^ 
1231 ]. In fact, an explanation of the spectrum of such states requires that the kernels used in formulating 
the associated bound-state problems are essentially nonperturbative, incorporating effects of DCSB 
which it has only recently become possible to express ^5; 36; ,37; 3^ . An equivalent formulation of the 
baryon bound-state problem is not yet available. Furthermore, technical difficulties associated with the 
analytic structure of rainbow-ladder kernels constructed using realistic interactions [l3; [HI; [11] have so 
far prevented computation of the spectrum of meson excited states, and the excited states and parity 
partners of ground-state baryons. 

Key elements in a successful spectrum computation are: symmetries and the pattern by which they 
are broken; the mass-scale associated with confinement and DCSB; and full knowledge of the physical 
content of bound-state kernels. These features are present in the informed use of a symmetry-preserving 
treatment of a vector x vector contact interaction. This underlies the success of that approach in Ref. 
which produced the first unified DSE description of the spectrum of light-quark meson and baryon 
ground- and excited-states, and is promising to provide a bridge between QCD and dynamical coupled- 
channels reaction models [8| . We undertake the current study in the expectation of similar reward in 
connection with strange hadrons. This is critical because contemporary hadron structure calculations, 
which typically omit meson-cloud effects, should not directly be compared with ex peri ment but instead 
with the bare- masses, -couplings, etc., determined via coupled-channels analvses [401: [4ll: [4^ [43| . 

In Sect. [3] we explain our interaction, and its application to mesons and colour- antitriplet quark- 
quark correlations. The latter bear no relation to the pointlike diquark degrees-of-freedom employed in 
some models of the constituent-quark type. Instead, they are dynamical elements that arise naturally in 

^ In connection with these phenomena, it is important to appreciate that the static potential measured in 
numerical simulations of quenched lattice-regularised QCD is not related in any known way to the question 
of light-quark confinement. It is a basic feature of QCD that fight-quark creation and annihilation effects are 
essentially nonperturbative and therefore it is impossible in principle to compute a potential between two light 
quarks. These points are elucidated, e.g., in Sect. 2.4 of Ref. [J]. 



Fig. 1 Poincare covariant Faddeev equation employed herein to calculate baryon properties. 'I' in Eq. (|29|l 
is the Faddeev amplitude for a baryon of total momentum P = pq + pd. It expresses the relative momentum 
correlation between the dressed-quark and -diquarks within the baryon. The shaded region demarcates the 
kern el of the Faddeev equation, Sect . 13.11 in which: the sing le line denotes the dressed-quark propagator, 
Sect. l2Al r is the diquark Bethe-Salpeter amplitude, Sect. l2.2l and the double line is the diquark propagator, 
Eqs. H33p . (|37p . Quarks within a diquark are correlated via gluon exchange but the kernel in this Faddeev 
equation expresses binding within the baryon through diquark breakup and reformation, which is mediated by 
exchange of a dressed-quark with momentum q. 

solving a Faddeev equation with QCD-based interactions; and, as we shall make clear, they are crucial in 
understanding the baryon spectrum. Section[3]describes the general structure of the Faddeev equations 
and solution amplitudes, and explains the impact of omitting resonant (meson cloud) contributions 
when constructing the Faddeev kernel. Our results are presented and discussed in Sect.Hl and we 
provide a summary and perspective in Sect. [5] 

2 Elements in the Faddeev Equation 

We base our description of baryon bound-states on a Poincare-covariant Faddeev equation, which 
is illustrated in Fig.[T] Introduced in Ref. its key elements are the dressed-quark and -diquark 
propagators, and the diquark Bethe-Salpeter amplitudes. All are completely determined once the quark- 
quark interaction kernel is specified and, as explained in the Introduction, we use 

g'^D^^ip -q) = Sf,^ ^^'^^^ , (1) 

where mc = 0.8 GeV is a gluon mass-scale typical of the one-loop renormalisation-group-improved 
interaction detailed in Ref. [22], and the fitted parameter am = 0.937r is commensurate with contem- 
porary estimates of the zero- momentum value of a running-coupling in QCD [i^; |4^; lU IH; IS US]- 
We embed Eq. ([T]) in a rainbow-ladder truncation of the DSEs, which is the leading-order in the most 
widely used, global-symmctry-prcserving truncation scheme |17{ . This means 

rAp,q)^i. (2) 

in the gap equation and in the subsequent construction of the Bethe-Salpeter kernels. 

One may view the interaction in Eq. ([Ij as being inspired by models of the Nambu-Jona-Lasinio 
type [m but our treatment is atypical. Used to build a rainbow-ladder truncation of the DSEs, Eqs. ([1]), 
^ produce results for low-momentum-transfer observables that are practically indistinguishable from 
those produced by more sophisticated interactions 0; H [HI; [i3l ■ 

2.1 Dressed-quark propagator 

Using Eqs. ([IJ, the gap equation for a quark of flavour / becomes 

5'^ (p) =i7-p + m/ + — — ^^7^S'/((7)7^, (3) 

where is the quark's current-mass. (Our Euclidean metric conventions are detailed in Add. El) 
Equation ^ possesses a quadratic divergence, even in the chiral limit. When the divergence is regu- 
larised in a Poincare covariant manner, the solution is 

Sfip)-^ = ij-p + Mf, (4) 
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Table 1 Computed dressed- quark properties, re quir ed a s inp ut for the Bethe-Salpeter and Faddeev equations, 
and computed values for in-hadron condensates [53 : [53 : 1541) . All results obtained with am ~ O.QStt and (in 
GeV) ylir ~ 0.24 , yluv = 0.905. N.B. These parameters take the values determined in the spectrum calculation 
of Ref. [il; and we assume isospin symmetry throughout. (All dimensioned quantities are listed in GeV.) 









Mo 


M„ 


Ms 


Ms/M^ 


1/3 


1/3 


1/3 


0.007 


0.17 


24.3 


0.36 


0.37 


0.53 


1.43 


0.241 


0.243 


0.246 



where M f is momentum-independent and determined by 



ds s 



s + Mf 



Our regularisation procedure follows Ref. [55|; i.e., we write 



AP 



dr c 



M2 



(5) 



(6) 



where rir.uv are, respectively, infrared and ultraviolet regulators. It is apparent from the rightmost 
expression in Eq. ([B]) that a finite value of rjr =: l/Ai,- implements confinement by ensuring the absence 
of quark production thresholds [l|; HI] . Since Eq. (|T]) does not define a renormalisable theory, then 
yluv ■= 1/tuv cannot be removed but instead plays a dynamical role, setting the scale of all dimensioned 
quantities. Using Eq. ([5]), the gap equation becomes 



Mf 



Mi 



4q;ir, 



(7) 



where C'^{(t)/(t = C (cr) = r{—l,aT^^) — r{ — l,aT?), with r{a,y) being the incomplete gamma- 
function. 

In Table [T] we report values of u- and s-quark properties, computed from Eq. ([7]), that will sub- 
sequently be used in bound-state calculations: the input ratio m^/m, where ffi = {mu + md)/2, is 
consistent with contemporary estimates [131 • N.B. It is a feature of Eq. ([7]) that in the chiral limit, 
ruf = mo = 0, a nonzero solution for Mq := lim,„^_j.oM/ is obtained so long as aiR exceeds a mini- 
mum value. With ylir,uv as specified in the Table, that value is af^ sa OAn. In the Table we also include 
chiral-limit and physical- mass values of the in-pseudoscalar-meson condensate [52l: [53t [S^. kh, which 
is the dynamically generated mass-scale that characterises DCSB. A growth with current-quark mass 
is anticipated in QCD [Tl: [ssl]. 



2.2 Mesons and diquark correlations 
2.8.1 Mesons 

The rainbow-ladder truncation of the gap and Bethe-Salpeter equations provides a good approximation 
for ground-state vector- and charged-pseudoscalar- mesons 0; IsH; HI; [6^ . We therefore employ it herein, 
in which case the homogeneous Bethe-Salpeter equation (BSE) for a meson comprised of quarks with 
fiavours /, 5 is 

rf,{k;P) = / 0^i,Sf{q + p)rr,{q;P)s,iqh,, (8) 

where P is the total momentum of the bound-state. This equation has a solution for P^ ~ ^^'fg' 
where nifg is the bound-state's mass. 

Here we illustrate the nature of the BSE via two relevant examples; viz., the negatively charged 
kaon and the kindred K* vector meson, both of which possess su flavour structure. The interaction in 
Eq. ((T)) supports a kaon Bethe-Salpeter amplitude of the form 



PxiP) - il5 Ek{P) + 771^757 -PFKiP) , 



(9) 
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wher(13 Mfi. ~ MsMu/[Ms + A/„]. If one inserts Eq. (|9]) into Eq. (|8l) and employs the symmetry- 
preserving regularisation of the contact interaction explained, e.g., in Ref. @, which requires 

,.1 



0= / da [C''iu;iMlMla,P^))+ Cl''{u:iMlMla,P^))] 



where 



uj{M^, Ml, a, P^) = Ml{\ - a) + aU] + a(l - a)P^ , 

CT{z) = ^z{dldz)C^{z) = z [r(0,MV2j - r(0, AfVj] , 

then the explicit form of the kaon BSE is 



EKiP) 

Fk{P) 



4aiR 



Aee -'yEF 



Ek{P) 
Fk{P) 



(10) 



(11) 
(12) 



(13) 



with 



q^K 



qy-K 

Afb 



q^K 

Aff 



C'-iioiMlMla^P^)) 



MuMs - a{\ - a)P^ - uj{MI, a, P^) 



c[^{u;iACMla,P'))y 



p2 



2M 



da 



R Ja 



{l-a)Mu + aMs 



C, {u:{Ml,Mla,P^)), 



2M| ^ 

p2 AbF; 



da 



MuMs + (1 - a)M;^ + aMl 



C, {u:{Ml,Mla,P^)). 



(14a) 
(14b) 
(14c) 
(14d) 



Equation is an eigenvalue problem, which has a solution for P^ = —mj^. The eigenvector 
is the kaon's Bethe-Salpeter amplitude, and in the computation of observables one must employ the 
canonically normalised amplitude; viz., the amplitude rescaled such that 

d 



1 = 



dP2 



nKiQ,P) 



2=P 



where 



7Tx(g,P) = 6trD 



(04^A'(-Q)^5.(<7 + P) PKiQ) 5„(<z) . 



Since Eq. ([T]) can only support a vector meson Bethc-Salpeter amplitude of the form 

Pk' =1^Ek*{P), 



where P^7^ = 0, the K* BSE is simpler; viz.. 



0, 



(15) 
(16) 
(17) 
(18) 



with 

-K 



2aiR 
37rm2, 



da 



MuMs - (1 - a)A/^ - - 2a(l - a)p2 



(c.(Af2,Af2,a,p2)), (19) 



where we have used Eq. (|10p . In this case the canonical normalisation condition can be written 



E 



1 _ 47raiR 

™2 ~ ^2 



(20) 



It should be plain that the analogous set of equations for the p-meson is obtained simply by replacing 
the s-quark by a d-quark throughout; and that for the by replacing the u-quark by a s-quark. Other 
states are discussed in App.IbI 



^ The choice one mcikes for the mass-dimensioned constant, Mr, has no effect on any result. 
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2.2.2 Diquark correlations 

The relevance of the rainbow- ladder meson BSE to the baryon Faddeev equation is explained, e.g., in 
Sect. 2.1 of Ref. 0; namely, in this truncation one may obtain the mass and Bethe-Salpeter amplitude 
for a colour-antitriplet quark-quark correlation (diquark) with spin-parity from the equation for a 
J~'^-meson in which the only change is a halving of the interaction strength [6l|. The flipping of the 
sign in parity occurs because it is opposite for fermions and antifermions. 

At this point it is appropriate to remark that the rainbow-ladder truncation generates asymptotic 
diquark states. Such states are not observed and their appearance is an artefact of the truncation. 
Higher-order terms in the quark-quark scattering kernel, whose analogue in the quark-antiquark chan- 
nel do not materially affect the properties of vector and flavour non-singlet pseudoscalar mesons, ensure 
that QCD's quark-quark scattering matrix does not exhibit singularities which correspond to asymp- 
totic diquark states [13; IH; E^l • Studies with kernels that exclude diquark bound states nevertheless 
support a physical interpretation of the masses, TO(gg)^p , obtained using the rainbow-ladder truncation; 
viz., the quantity ^(^ggyj^ ^/''^(qq)jp n^&Y be interpreted as a range over which the diquark correlation 
can propagate before fragmentation. 

This caveat expressed, one may write the contact-interaction rainbow-ladder BSE for a colour- 
antitriplet diquark constituted from quarks with flavour /, g: 



(21) 



where C is the charge-conjugation matrix, Eq. (|A.9p . 

Capitalising further on the connection between the meson and diquark sectors, one may readily 
write explicit forms of the BSEs and canonical normalisation conditions for scalar ([/g]) and axial- 
vector ({//}, {fg}) diquark correlations. For example, the Bethe-Salpeter amplitude for a = 0+ 
[s, u]-diquark is 



^[iu(^)=^75i?M„,(P) + 
which satisfies the following BSE 



1 



l5l-PF[su]„^ (P), 



2qir 



Ab_e Aef 

qcK q^K 



In this case the canonical normalisation condition is 



(22) 



(23) 



(24) 



where 



n^suu {Q, P) = 4trD j -0^,r[suu i-Q)^^sM + p) Pisu^ (Q) Su 



(25) 



Compared with Eq. ([T5|) , the colour factor is different owing to the fact that diquarks are colour- 
antitriplets not singlets. 

Following this pattern one may immediately write the BSE for — 1^ {su} diquark correlations; 
viz., 

^-lK^'{-mlu}) = 0, (26) 
and the canonical normalisation condition: 

1 



{iStl} 



(27) 



The analogous set of equations for axial-vector {uu}- and {wd}-diquarks arc obtained by replacing 
the s-quark by either a u- or d-quark throughout; and that for the {ss}-diquark by replacing the 
u-quark by a s-quark. Other correlations are discussed in App.lB] 
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Table 2 Row 1: Quark-core masses of ground-state mesons computed using our symmetry-preserving regular- 
isation of the vector x ve ctor contact interaction, with the input from Table[T| Row 2: Except for scalar mesons, 
values drawn from Ref. [G^, with weighted averages of mass-squared values reported, where appropriate. For 
the isoscalar-scalar meson we list an estimate for the state's dressed-quark core [63; [gJ . Nothing is known 
about this value for the 1 = 1/2 scalar. Rows 3 and 4 repeat this pattern for the mesons' first radial excitation. 
The theory error in Row 3 displays the outcome of varying the location of the node in the radial excitation's 
Bethe-Salpeter amplitude: 1/dj — 2M^{1 ± 0.2). An asterisk-marked mass in Row 4 indicates a state whose 
properties are poorly determined. (All dimensioned quantities are listed in GeV.) 







niK 


mp 












rriKi 




n=0 DSE 


0.14 


0.50 


0.93 


1.03 


1.13 


1.29 


1.40 


1.38 


1.48 


1.59 


expt. 


0.14 


0.50 


0.78 


0.89 


1.02 


1.0 - 1.2 




1.23 


1.34 


1.42 


11=1 DSE 


1.33±o.o6 


1.33±o.o7 


1.29±o.o5 


l-40±o.o5 


1-51±0.05 


1.42±o.02 


1.53±o.o2 


1.47±o.o2 


l-57±o.oi 


l-67±o.o2 


expt. 


1.3±o.i 


1.46* 


l-46±o.o3 


1.68* 


1.68±o.o2 






l-65±o.o2 







2.2.3 Mesons: computed masses 

Before discussing the results presented in Table [H it is necessary to recapitulate on an important 
modification of the rainbovif-ladder Bethe-Salpeter kernel that one should implement before prediction 
and comparison with experiment are meaningful. It has long been known that the rainbow-ladder 
truncation describes vector meson and flavour-nonsinglet pseudoscalar-meson ground-states very well 
but fails for their parity partners j2^ : [23l : l65l : i66; 67; 6g|. The origin and solution of this longstanding 
puzzle are now available following a novel reformulation of the BSE [s^ , which is valid and tractable 
when the quark-gluon vertex is fully dressed. In employing this approach to study the meson spectrum 
it was found that DCSB generates a large dressed-quark anomalous chromomagnetic moment and 
consequently that spin-orbit splitting between ground-state mesons is dramatically enhanced [36t [btI : 
[69| . This is the mechanism responsible for a magnified splitting between parity partners; namely, there 
are essentially nonperturbative DCSB corrections to the rainbow-ladder kernels, which largely-cancel 
in the pscudoscalar and vector channels but add constructively in the scalar and axial- vector channels. 

With this in mind, we follow Ref. Q and introduce spin-orbit repulsion into the scalar- and 
pseudovector-meson channels through the artifice of a phenomenological coupling g^Q < 1, introduced 
as a single, common factor multiplying the kernels defined in Eqs. (jB.Qp . ()B.15p . The valu^ 

5SO = 0.24 (28) 

is chosen so as to obtain the experimental value for the ai-p mass-splitting, which we know to be 
achieved by the corrections described above [s^ [s^ [6^. It is noteworthy that the shift in niai is 
accompanied by an increase of ma- and that the new value matches an estimate for the gg-component 
of the cr-meson obtained using unitarised chiral perturbation theory (63l : |63 | . 

This expedient produces the results for scalar and axial- vector mesons in Rows 1,3 of Table [2j which 
reports calculated results for meson masses and compares them with available empirical values. This 
information is represented pictorially in the left panel of Fig. [5] Owing to our choice for the current- 
quark masses, to^ and niK agree with experiment. All other computed values for ground-states are 
greater than the empirical masses, where they are known. This is typical of DCSB-corrected kernels that 
nevertheless omit resonant contributions; i.e., do not contain effects that may phenomenologically be 
associated with a meson cloud. In Table[3]we list the canonically normalised Bethe-Salpeter amplitude 
for each meson. These are the quantities used in calculating observable properties of mesons; and 
comparison with the kindred diquark amplitudes, listed in Table [5l will subsequently be instructive. 

In practical calculations, meson-cloud effects divide into two distinct types. The first is within the 
gap equation, where pscudoscalar meson loop corrections to the dressed-quark-gluon vertex act to 
reduce uniformly the mass- function of a dressed-quark [ll; [t^; [Zll; III; [ill ■ This effect can be pictured 
as a single quark emitting and reabsorbing a pseudoscalar meson. It can be mocked-up by simply 
choosing the parameters in the gap equation's kernel so as to obtain a dressed-quark mass-function 



^ NB. gso ~ 1 means no spin-orbit repulsion. The mass changes slowly with diminishing gso; e.g., gso = 0.50 
yields = 1.23 GeV. 
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n K p K* 



1.2 

> 1. 

CD. 0.8 



C 



0.6 



tO.4 
" 0.2 



1 



7T K p (p a- K 



Fig. 2 Left panel: Pictorial representation of Table [2] Circles - computed ground-state masses; squares - 
computed masses of radial excitations; diamonds - empirical ground-state masses in Row 2; and triangles - 
empirical radial excitation masses in Row 4. Right panel: Circles - computed splitting between the first radial 
excitation and ground state in each channel; and triangles - empirical splittings, where they are known. The 
dashed line marks a splitting of 0.1 GeV. 



that is characterised by a mass-scale of approximately 400 MeV. Such an approach has implicitly been 
widely employed with phenomenological success S E [Zl III ■ We employ it herein. 

The second type of correction arises in connection with bound-states and may be likened to adding 
pseudoscalar meson exchange between drcssed-quarks within the bound-state [t^ [t^; ^Bh IzS HO; HH El] ; 
as opposed to the first type of effect; i.e., emission and absorption of a meson by the same quark. 
The typc-2 contribution is that computed in typical evaluations of meson-loop corrections to hadron 
observables based on a point-hadron Lagrangian. These are the corrections that should be added to 
the calculated results in Tabled The most complete computation of this sort predicts that such effects 
reduce rUp by 0.13 GeV [U. Applied to our result, this would produce mi^oop-coTected ^ q.SGcV, in 
good agreement with the empirical value of 0.78 GeV. 

These observations underpin a view that bound-state kernels which omit type-2 meson-cloud cor- 
rections should produce dressed-quark-corc masses for hadron ground-states that are larger than the 
empirical values. As we shall see, this is uniformly true herein. Moreover, thisperspective also has 
implications for the description of elastic and transition form factors d; [2^; HI; Hjl . 

The situation for radially excited states is less clear. This may be seen from the right panel of 
Fig.[2j which depicts the computed mass splitting between ground-states and the first radial excitation 
in each channel; and also provides a comparison with experiment, when that is available. The compar- 
ison suggests that our formulation of the contact interaction kernels produces the correct trend but 
underestimates the splitting by ^ 0.2 GeV. We note that this mismatch is reduced if type-2 meson- 
cloud corrections to the masses of radial excitations are smaller than for ground-states. On the other 
hand, it might simply be that this underestimate is an error arising from the expedient we employ 
in order to define radial excitations within the contact interaction framework, which is discussed in 



Table 3 The structure of meson Bethe-Salpeter amplitudes is described in Sect . 12. 2!T] and App.jB] Here we list 
the canonically normalised amplitude associated with each of the BSE eigenstates in Table[2] Only pseudoscalar 
mesons involve two independent amplitudes when a vector x vector contact interaction is treated systematically 
in rainbow-ladder truncation. 







TUk 


nip 


niK' 




ma- 


771k 


777ai 


771 K J 




n=0 


Egg 

Eqq 


3.60 
0.48 


3.86 
0.60 


1.53 


1.62 


1.74 


OA? 


0.47 


0.31 


0.31 


0.31 


n=l 


Eqq 

Cqq 


0.83 
0.05 


0.76 
1.18 


0.72 


0.70 


0.66 


0.34 


0.35 


0.28 


0.28 


0.28 
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S # # 

♦ 

♦ • 
♦ • 

♦ 



[udjQ. [usio* luulr lusir Issl,. ludio- luslo- luuir lusl,- {ssl,- 

n K p K* ^ a- K f-^ 




« 0.25 



[udlo* [uslo* {uuir {ush* {SSI,. ludIo- luslo- luuIr lusIr Issl,- 



n K p K* 



cr K a-i K-i 



Fig. 3 Left panel: Pictorial representation of Table [l] Diamonds - ground-state diquark masses in Row 1; 
circles - ground-state meson masses in Row 2; triangles - masses of diquark first radial excitations in Row 3; 
and squares - masses of meson radial excitations in Row 4. Right panel: Diamonds - for diquarks, computed 
splittings between first radial excitation and ground state; and circles - for mesons, computed splitting between 
the first radial excitation and ground state in each channel. The dashed line marks a splitting of 0.1 GeV. 



Add. IB.i1 Given that possibility, one must allow that our predictions for the dressed-quark-core masses 
of hadron first radial excitations might be ~ 0.2 GeV too small. 



2.2.4 Diquarks: computed masses 

The preceding discussion of systematic trends within our predictions for meson masses is important 
to understanding the results of our Faddeev equation studies because of the connection between the 
meson and diquark BSEs, outlined in Sect. 12. 2T^ predictions for mesons masses determine the diquark 
spectrum and hence impact heavily on the baryon spectrum. In Table S] we therefore present results 
for the masses of diquark correlations and compare them with the meson masses in Table [D This 
information is also depicted in the left panel of Fig.[3l 

It is plain from Fig. [3] that the level ordering of diquark correlations is precisely the same as that for 
mesons. Moreover, in all diquark channels, except the scalar, the mass of the diquark's partner meson 
is a fair guide to the diquark's mass: the meson mass bounds the diquark's mass from below; and the 
splitting is always less than 0.13 GeV and decreases with increasing meson mass. 

The scalar diquark channels are particular owing to DCSB and the Goldstone boson character of 
the partner pseudoscalar mesons. We note that in a two-color version of QCD, the scalar diquark is 
also a Goldstone mode [s^ [s^, a long-known result of Pauli-Giirscy symmetry [s^; H^. A memory 
of the symmetry persists in the three-color theory and is evident here in low masses for the scalar 
diquarks. That they are low is highlighted by the right panel of Fig. [31 which shows large splittings 



Table 4 Row 1: Quark-core masses of diquark correlations that play a role in the octet and decuplet spectra of 
baryons, computed using our symmetry-preserving regularisation of the vectorx vector contact interaction, with 
the input from Table[T] Row 2: Ground-state meson masses from Row 1 of Table[21 Rows 3, 4 repeat the pattern 
of Rows 1, 2 for the diquarks' first radial excitations. The theory error in these rows displays the outcome of 
varying the location of the node in the radial excitation's Bethe-Salpeter amplitude: 1/djr = 2M^(1±0.2). (All 
dimensioned quantities are listed in GeV.) 





[u, d\o+ 


[s,u\o+ 


|u,w|i+ 




|s,s|i+ 


[u, d\o- 




{u,u\i- 






n=0 


m 


0.78 


0.93 


1.06 


1.16 


1.26 


1.37 


1.47 


1.45 


1.55 


1.65 




qq 


0.14 


0.50 


0.93 


1.03 


1.13 


1.29 


1.40 


1.38 


1.48 


1.59 


n=l 


qq 


1.34±o.o5 


1.35±o.o5 


1.32±o.o4 


1.42±o.04 


1.53±o.o4 


1.48±o.o3 


1.57±o.o2 


1.52±o.oi 


1.62±o.o2 


1.71±o.ui 




qq 


1.33±o.o6 


1.33±o.o7 


1.29±o.05 


1.40±o.o5 


1.51±0.05 


1.42±o.o2 


1.53±o.o2 


1.47±o.o2 


1.57±o.oi 


1.67±o.o2 
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Table 5 The structure of diquark Bethe-Salpeter amplitudes is described in Sect. 12. 2!2] and App.[Bl Here we 
list all canonically normalised amplitudes that are relevant to the baryons we consider. Only scalar diquarks 
involve two independent amplitudes. 





[u, dJo+ 


ls,w|o+ 






|s,s|i+ |u,dJo- 




{u,u\^- 




{■s,s|i- 


Eqq 


2.74 


2.91 


1.30 


1.36 


1.42 0.40 


0.39 


0.27 


0.27 


0.26 


Fqq 


0.31 


0.40 

















between the ground and excited states in the scalar diquark channel. Notwithstanding this, the scalar 
diquark correlations are split widely from the true Goldstone mode mesons. 

In constructing baryon Faddeev equations, the canonically normalised diquark Bethe-Salpeter am- 
plitudes are critical because they determine the strength of the correlation within a given baryon. We 
list them in Table [S] Notably, the amplitudes of the positive-parity states are much larger than those 
of states with negative-parity: for like-flavour content, the ratio always exceeds 5. This pattern repeats 
that established by the ground-state mesons, see Table [31 



3 Baryon Faddeev Equations 

3.1 General structure of the Faddeev amplitudes 

A spin- 1/2 baryon is represented by a Faddeev amplitude [i^] 

•Z/ = 1^-1 + IP's + "^3 , (29) 

where the subscript identifies the bystander quark and, e.g., ^1^2 are obtained from 1^3 by a cyclic 
permutation of all the quark labels. We employ the simplest realistic representation of •P', so that an 
octet baryon is composed from a sum of scalar and axial- vector diquark correlations: 

^:,{p,,a„^,) = ^t,+9tt, (30) 

with {pj,aj,ipj) the momentum, spin and flavour labels of the quarks constituting the bound state, 
and P = pi + P2 + P3 the system's total momentum. 

It is conceivable that pseudoscalar and vector diquarks could play a role in the Faddeev amplitudes 
of ground-state = (1/2)"*" baryons. However, such correlations have opposite parity and hence can 
only appear in concert with nonzero quark angular momentum. Since one expects ground-states to 
possess the minimum possible amount of quark orbital angular momentum and these diquark corre- 
lations are significantly more massive than the scalar and axial- vector (see Fig.lSj, they can safely be 
ignored in computing properties of the ground state. 

In order to assist in explicating the structure of the diquark pieces in Eq. pop , we define a set of 
flavour matrices 





1 




1 







^l = [ud] _ 


-10 


^2=[us] _ 


00 


.(.3=[tis] _ 


1 




00 




-10 




0-10 




V2 0' 




"0 1 0' 






"0 r 


^4={uu} 


00 


^ ^5={ud} ^ 


1 


.(.6={ms} _ 









00 











1 




"0 0' 




"0 0' 






00 1 ' 


^7={dd} _ 


y/2 


.(.8={!is} _ 


1 


^9={ss} ^ 




00 









1 






OOy/2 



Employing these matrices, the scalar diquark piece in Eq. pop can be written 



(32) 
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where: K =pi+p2 ^:p{i2}, P[i2] = Pi - P2, t (-P{i2} + 2p3)/3; 

1 



2I0+ 



(33) 

is a propagator for the scalar diquark formed from quarks 1 and 2, with tiii^h^j]^^ the mass-scale 
associated with this [(^i(/?2]o+ diquark; ^j^^^^^j is the canonically-normalised Bcthe-Salpeter amplitude 
describing the relative momentum correlation between the quarks; 5, a 4 x 4 Dirac matrix, describes 
the relative quark-diquark momentum correlation within the baryon; and the spinor satisfies 



(34) 



with Mxp the baryon mass obtained by solving the Faddeev equation. Wc note that vF also possesses 
another column- vector degree of freedom; viz., 





[ud]u 




[us]u 




[us]s 




Up = 


{uu}d 


, us+ = 


{uu}s 




{us}s 






{ud\u 




{us}u 




{ss}u 



y/2[ud]s 
[ud]s — [ds]u 
{us}d — {ds}u 



(35) 



Owing to our assumption of isospin symmetry, the unlisted octet charge states are degenerate with 
their listed partners. 

The axial-vector part of Eq. ((30)) is 



where 



1 



{VlV2}i+ 



(36) 
(37) 



is a propagator for the axial-vector diquark formed from quarks 1 and 2 and the other elements in 
Eq. pop are obvious analogues of those in Eq. ([5^ . 

In connection with decuplct baryons we note that it is not possible to combine an isospin-zero 
diquark with an isospin-1/2 quark to obtain isospin-3/2 and hence the A is comprised solely from 
axial- vector correlations. This sets the pattern for the remaining dccuplet baryons, which may therefore 
be expressed via 

(38) 



with 
2)^K,(pj,aj,^j) 



(39) 

where (P) is a Rarita-Schwingcr spinor and, as with octet baryons, in constructing the Faddeev 
equations we focus on that member of each isospin multiplct which has maximum electric charge; viz.. 



{uu}s 
{us}u 



{us}s 
{ss}u 



un = [{ss}s] 



(40) 



The general forms of the matrices 5*(^;P), S^{t]P) and 'D^p'{i;P), which describe the 
momentum-space correlation between the quark and diquark in the octet and decuplet baryons, re- 
spectively, are described in Refs. [H; |8^. The requirement that S'^{l;P) represent a positive energy 
baryon entails 

S'^{e;P) = st{e■,P)Iu+[^l■i-i■PIB) st{t,P), (41) 
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where (iD)rs = Srs, ~ 1, = —1. In the baryon rest frame, describe, respectively, the upper, 
lower component of the bound-state baryon's spinor. Placing the same constraint on the axial-vector 
component, one has 

6 

J^{e;P) = J2 Vl P) 75 AZ {t,P), (42) 

where = ^ i • P P^, 7^ = 7^ -f 7 • P P^) 

Al=^-i^P,, Al = -iP,, Al^^-i^i^, 

At - Al^ 7.^ -Al, Al^ ■ - At . ^^'^ 

Finally, requiring also that 2?^^ {£; P) be an eigenfunction of A^{P), one obtains 

(^; P) = ^) + 75<" (^; p) < , (44) 

with 5"^ and given by obvious analogues of Eqs. (j4ip and (|42|) . respectively. 

Having in hand detailed forms for the dressed-quark propagators, the diquark Bethe-Salpeter am- 
plitudes and the diquark propagators, it is now possible to write Faddeev equations for the baryons. As 
apparent in Fig.[TJ the kernels of those equations involve diquark breakup and reformation via exchange 
of a dressed-quark. In proceeding we follow Ref. @ and make a drastic simplification; namely, in the 
Faddeev equation for a baryon of type P, the quark exchanged between the diquarks is represented as 

S^{k) ^ ^ , (45) 

where the superscript "T" indicates matrix transpose, / is the quark's flavour and qb is discussed below. 
This is a variant of the so-called "static approximation," which itself was introduced in Ref. [l^l and 
has subsequently been used in studies of a range of nucleon properties @; |9l| . In combination with di- 
quark correlations generated by Eq. ([1]), whose Bethe-Salpeter amplitudes are momentum-independent, 
Eq. (j45[) generates Faddeev equation kernels which themselves are momentum-independent. It follows 
that Eqs. (|^T|) . simplify dramatically, with only those terms that are independent of the relative 
momentum surviving: 

5*(£; P) ^ 5*(P) = .s*-(P) Id , (46a) 
P) ^ SlliP) = af (P) Z757M + 4(^^)75^^ • (46b) 



3.2 Explicit example: A baryon 

In App. lC.ll we derive the Faddeev equation for the resonance, whose simple spin-flavour structure 
makes it a useful illustrative example; and in App. lC.2l the nucleon's Faddeev equation is reported, with 
indications of the minor changes in analysis that are necessary to complete its derivation. Here, on the 
other hand, we illustrate different aspects of the construction of the Faddeev equation by considering 
the A baryon, which is an isospin-zero, J = (1/2)+ state constituted from a single quark of each flavour 
and hence has a complicated spin- flavour amplitude. 

Recalling comments made when opening Sect. 13. ll five possible diquark combinations arc possible 
for the ground-state A: 

s[ud]a+, d[Ms]o+, w[(is]o+, d{tts}i+, u{ds}i+. (47) 

Of these, s[wd]o+ has / = 0, whilst the others do not possess good isospin. This lies behind a mixing that 
obscures the A and octet baryon isospin-eigenstates and entails that building the flavour structure 
of the A Faddeev kernel is quite involved. States of good isospin can be constructed as follows: with 



/ [ud]s \ 




/I 














\ 


[us\d 







1 

V2 


1 

72 










[ds\u 


,o = 





1 

72 


1 

\/2 










{us}d 













1 

x/2 


1 

x/2 




\{ds}u) 












1 

V2 


1 

\/2 


/ 
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then each of the entries in the new column vector 



V 



1 ^J2[ud]s, 


I = 




[us\d — [ds\u , 


I = 





[us\d + [ds\u , 


I = 


1 


{us}d — {ds\u 


,1 = 





\ {us}d + {ds\u 


,1 = 


1/ 



(49) 



is an eigenvectors of isospin, with the eigenvalue indicated. 

A consideration of Fig.[T] reveals that the column vector V satisfies a Faddccv equation of the form 
V = KudsV , which it is helpful to write explicitly: 



[us]d 
[ds]u 
{us}d 
\{ds}u/ 



I 





us\ , [ud] 



K 

\us\ 



i^{ud\\us\ i^\ud\\ds\ i^\u.d\,{us} i^\ud\\ds} \ 



us] , [ds\ 





K 



{us} , [ud\ 



K 





[ds] ,{i/.s} 



us] ,{c/s} 





\K{ds}.\ud] K{ds}, 



{t(s},[ds] 









K 



K. 



{(is},{us} 



{ us} ,{ds} 





[us\d 
[ds\u 
{us\d 
\{ds}uj 



(50) 



where, e.g., -ftr[„d].[„s] describes the breakup of a [us] scalar diquark through emission of a dressed u- 
quark, which then joins the d-quark to form a [ud] scalar diquark, leaving the s-quark as a bystander. 
In the kernel, the repeated flavour label always indicates the exchanged quark. Using the notation 
connected with Eq. pip , one may write 



Wc now make two entries explicit: 



/ Ki2 Kie Kis\ 
K21 K23 K28 
K31 K32 Kse 
Kei Ke3 K^s 

\Ks1Ks2 Kse J 



(51) 



K12 = Upt r[us]ilds)Si^t r[ud]i~~kud)Sdild)'^[us]{lus)UF ■ 
= r[us]ilds)S'^r[ud]i-kud)Sd{ld)^[us]{lus) =■ ^2, 



Kl3 = Ujt^r[ds]ilds)S^t'-'- r[ud]i-kqq)Suilu)'^[ds]ilds)UF, 



iT^lT 



-r[ds]{lds)Sjr[ud]i-kqq)Suilu)^[ds]ilds 



-%-3 



(52a) 
(52b) 
(52c) 
(52d) 



where Uj = (1, 1, 1) is a device that we use to collapse the flavour structure. Repeating this procedure 
for all entries and using isospin symmetry, one arrives at 



^uds 



/ ^12 —"KlS ^16 \ 

'K2I 0(^23 X,28 

-%.2l 'K.23 ^^28 

~^61 %_io2 ^68 

\ %Jo\ 'K&2 %ja8 / 



(53) 



As the kernel K^ds contains mixing that obscures the A and Z"*^ isospin-eigcnstate baryons, wc now 
employ the matrix O in Eq. (|48p and construct the non-mixing kernel: 



%,ds^OK,,dsO' 



1 







- 




\ 




V2f^21 


— 'K.23 





-K28 













'K.23 





K28 






-'K&3 





-^68 





\ 








^63 





^68 / 



(54) 
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Rows 1, 2, 4 map 1 = into itself, whereas rows 3, 5 do the same for {I,Iz) = (I7O). Focusing on 
the / = sector, one arrives at the following Faddeev equation for the A baryon: Va = ^Ws^i i-^-' 
explicitly. 



Va 



1 

7^ 



/ y/2[ud]s 

[us]d — [ds]u 
{us}d — {ds}u 



72^^12 -V2i^i6 




/ 72 [ltd] s 

[us]d — [ds]u 
y {us}d — {ds}u 



(55) 



Acknowledging that the axial-vector diquark components of all baryon amplitudes involve two 
scalar functions, as evident in Eq. (j46bp . we re-express as follows: 



/ 









-\/2^?Cf62 






V2^21 


n'A 
-A.23 


Ht'A 
-A.281 


n^A 
A.282 






-a/2^^6i1 


n'A 
A.6i3 


n'A 
A.6181 


n^A 
A.6182 




V 


-V2^ii 


qy-A 
A.623 


n(-A 
A.6281 


A.6282 





(56) 



i.e., adding subscripts associated with the diquark labels 6,8. Having arrived at this form, one can 
write compact algebraic expressions for each of the entries, which we list in Adp.[U31 

Equation ([55)) is an eigenvalue problem whose solution yields the mass for the dressed-quark-core 
of the yl- resonance. If one sets gyi = 1, then the Faddeev equation yields ruA = 1.39 GeV and the 
unit-normalised eigenvector 



[6,8] 
A2 

0.65 0.59 -0.47 -0.020 



„[2,3] [6,8] 



a 



(57) 



The eigenvector corresponds to the following probabilities: P\ = 43%, pjj^''^' = 35%, pj^^'*' = 22%; 



i.e., the ground-state A is 22% axial vector diquark and 78% scalar, with a roughly equal probability 
of both / = scalar diquark configurations. Comparison with Eq. (jC.18p reveals an interesting result; 



P 



[6,8] 



p4 



p5 



So, despite the differences in flavour structure, the net axial- vector diquark 



content of the nucleon and A is the same. 



3.3 Pion-loops and baryon masses 

We discussed the nature of resonance contributions to meson bound-state kernels on page [T] Similar 
remarks may be made in connection with baryon Faddeev equations and this provides us with a means 
by which to fix qb in Eq. ps)) : viz., formulae such as those in Ref. [92| can be used to estimate the size 
of typc-2 meson-loop corrections to baryon masses computed using the Faddeev equations herein. The 
straightforward application of such expressions, using a common mcson-baryon form-factor mass-scale 
of 0.8 GeV, yields a shift of (—300 MeV) in tun and (—270 JVIeV) in niAi from which one may infer that 
our type-2 uncorrected Faddeev equations should produce rriN = 1.24 GeV and ruA = 1.50 GeV, values 
which Adds. [U^IIC. 21 show to be of the appropriate size. For the Z\-resonance there is another estimate, 
which is arguably more sophisticated. Namely, that produced by the Excited Baryon Analysis Center 
(4^ . which used a realistic coupled-channels model to remove meson dressing from the A and expose 
a dressed-quark-core bare-mass of 1.39 GeV. Following these observations, which are condensed from 
Ref. 0, we return to Eq. PSI) and choose 

38 gB=N,A,E,E = 1-18 , .910 := gA,i:-,S'j2 = 1-56 , (58) 

so that the Faddeev equations in Adds. [U~1IC.2I produce m^v = 1.14 GeV, niA = 1.39 GeV, Sm = 
niA — rriN = 0.25 GeV because these outcomes are consistent with the information presented above 
and Refs. [S^H^illJ- N.B. Notwithstanding this common choice for the couplings within the octet and 
separately within the decuplet, the kernels in different strangeness sectors are still distinct owing to 
the appearance of the dressed-quark mass in the denominator of Eq. (j45p . 
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Table 6 Dressed-quark-core masses of ground-state octet and decuplet baryons, those of their radial excita- 
tions and of all their parity partners, comp uted using our symmetry- preserving formulation of the vector x vector 
contact interaction supplemented by Eqs. (|58|) and the material in Sect. [21 Row 1: Baryon ground-states. The 
lowest mass state in each channel has positive parity. Row 3: First radial excitations of the ground-states, 
which provide the second level in each channel. The theory error in this row displays the outcome of varying 
the location of the node in the radial excitation's Faddeev amplitude: l/dj- = 2M^(1 ± 0.2). Row 5: Parity 
partners of the baryon ground-states, which provide the third level in each channel. Row 7: First radial excita- 
tions of the parity partner to each of the baryon ground-states, which provide the fourth level in each channel. 
Masses in the rows labelled "expt." are taken from Ref. [SJ. To explain the notation we observe that the entry 
1.44p^^ in the "N" column is read as "A''(1440)Pii" , and a hyphen in any position indicates that no empirically 
known resonance can confidently be associated with the associated state. (All dimensioned quantities are listed 
in GeV.) 
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DSE 


1.14 


1.26 


1.35 


1.51 
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1.39Pi3 


1.53pi3 
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DSE 
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1.44p,, 


1.89±o.o7 
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2.09±o.oi 
1.66Pii 
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I.6OP33 
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2.04±o.o4 


2.14±o.04 
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DSE 


2.30 


2.40 


2.41 
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2.33 


2.44 


2.55 
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expt. 


1.54sii 


1.67soi 


1.75sii 




1-701533 


1.67di3 


1.82d^3 




p = 




a— 


1 


DSE 
expt 


2.35±().oi 
1.65sii 


2.45±o.oi 
1.80soi 


2.46±o.oi 


2.55±o.oi 


2.39±o.oi 


2.49±o.oi 
1.94ci3 


2.59±(,.oi 


2.70±o.oi 



4 Baryon spectrum 

4.1 Faddeev equations 

The Faddeev equations for all ground-state baryons in our formulation of the contact interaction are 
presented in App.lCl 

Turning to radial excitations we note that, in analogy with mesons, the leading Tchebychev mo- 
ment of the bound-state amplitude for a baryon's first radial excitation should possess a single zero. 
Hence, as advocated in Ref. @, it is possible to estimate masses for these states by employing the 
expedient described in App. IB.l1 With the zero located as prescribed in Eq. (jB.ip . no new parameters 
are introduced. Following this reasoning, the Faddeev equation for the first radial excitation of each 
baryon is simply obtained by making the following replacement throughout the equations in App.lCl 

Ci(<7)^7i((t). (59) 

In a more general setting one might imagine that a baryon's first radial excitation could be an 
admixture of two components: one with a zero in the Faddeev amplitude, describing a radial excitation 
of the quark-diquark system; and the other with a zero in the diquark's Bethe-Salpeter amplitude, 
which represents an internal excitation of the diquark. The procedure in App. IB.H can conceivably 
distinguish between these components via a mixing term whose strength is cx ^(/c()^j5£'(/g)'p 1 where 
the latter is the excited diquark's amplitude. Owing to orthogonality of the two-body ground- and 
first-radially-excited states, we anticipate that this mixing term is negligible. Under this assumption, 
a baryon's first radial excitation is predominantly a radial excitation of the quark-diquark system. 
Should a state constituted from a radially-excited diquark exist, then its mass will be larger because 

^ U9)jp' (fg)jp' ^ (fg)jp' {fg}jp' 

The case of parity partners is more complicated. One must typically reanalyse each isospin/flavour 
channel, beginning with the equation represented by Fig.[Tl because of the altered Dirac structure 
of the diquark correlations that arc involved. For example, the nucleon's parity partner is composed 
of pscudoscalar and vector diquark correlations and its Faddeev amplitude must change sign under 
a parity transformation. These alterations lead to changes in the locations of the 75 matrices in the 
Faddeev equation and also a reduction in the number of terms because the pscudoscalar diquark 
does not possess a i^(P)-component [Eq. (j^^ cf. Eq. (|B.13|) ]. With the exception of the Z\(| ), whose 
Faddeev equation is written in Eq. (IC.12p . we do not report parity-partner Faddeev equations herein. 



16 




N1N3N5N7 A1A3A5A7 BiBgBsBy A1A3A5A7 3;3*3*h^ n^n^^j 



Fig. 4 Left panel: Pictorial representation of octet masses in Table [S] Circles - computed masses; and 
diamonds - empirical masses. On the horizontal axis we list a particle name with a subscript that indicates 
its row in the table; e.g., A'^i means nucleon column, row 1. In this way the labels step through ground-state, 
radial excitation, parity partner, parity partner's radial excitation. Right panel: Analogous plot for the decuplet 
masses in Table [6l 



As will be seen from Eq. (jC.12p . however, the Faddeev kernels that one obtains via this reanalysis have 
similar forms to those for the positive parity states, listed in App.jCj but sign changes are introducecfl 
and, throughout, one has the replacements 

^(/9),p E{fa)j-p ' "^(/9),,p ^ "^(/9),,-p- (60) 
Equations for the radial excitations of the parity partners are readily obtained by using Eq. ([5^ . 



4.2 Computed Results 

We have solved all the Faddeev equations, supplemented by Eqs. ([55]) and the material in Sect. [21 and 
obtained the masses and eigenvectors of the ground-state octet and decuplet baryons. The masses are 
listed in Table [6l with this information also depicted in Fig.|4l 

It is evident in Fig. 2] that, as with mesons, our computed baryon masses lie uniformly above the 
empirical values. We view this as a success because our results are those for the baryons' dressed-quark- 
cores, whereas the empirical values include effects associated with type-2 meson-cloud effects, which 
typically produce sizable reductions [40l : |42| . Our values may reasonably be viewed as bare-mass inputs 
appropriate for dynamical coupled-channels analyses of the hadron spectrum Q. This was explained 
and illustrated for the nucleon and Zi-resonance in Sect. 4.5 of Ref. [1] and in particular for the Roper 
resonance in Ref. Q . Here we reiterate those instances in which a comparison can be made: 
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^1232^33 


^1700^33 


Table El (DSE) 


1.14 
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2.30 


2.35o.oi 


1.39 


2.33 


M% (Ref. [40j) 


1.24 




2.05 


1.92 


1.46 


2.25 


M% (Ref. [42j) 




1.76 


1.80 


1.88 


1.39 


1.98 



(61) 



where M^, when it appears, is the relevant bare mass inferred in the associated coupled-channels 
analysis. These bare masses have hitherto been uncertain and dependent on model details. However, 
as wc made no attempt to fit them, their proximity to our results suggests that it might now be possible 
to place these bare masses on firmer ground, investing them with meaning within the context of hadron 
structure calculations that have a traceable connection with QCD. 

In Fig.[S]wc plot the theory-experiment mass differences computed, where possible, from Table [51 
This difference is uniformly less-than 0.2 GeV for ground-states. Moreover, for decuplet ground-states 
it decreases uniformly with increasing strangeness, a result that owes to the simplicity of the decuplet 

* These changes were inadvertently overlooked in Ref. 0|. They lead to increases in the mass of parity 
partners but the effect is always smaller than 4%. 
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Fig. 5 Theory-experiment mass differences computed, where possible, from Table [6l The dashed line is a 
difference of 0.2 GeV; the dotted line, 0.37 GeV; and the dot-dashed line, 0.68 GeV. Horizontal axis: particle 
name with a subscript that indicates its row in the table; e.g., A''i means nucleon column, row 1. In this way 
the labels step through ground-state, radial excitation, parity partner, parity partner's radial excitation. 



Faddeev amplitudes. For octet radial excitations the difference is roughly 0.4 GeV; and for higher ex- 
citations it is around 0.7 GeV. This identification of band clustering for the excited states is. perhaps, 
arguable but the scatter can be attributed, at least in part, to experimental uncertainty. Notwithstand- 
ing this, these features suggest a regularity in the strength of typc-2 meson-loop corrections within 
excitation bands, with the case being strongest for ground-states. 

In this connection we remark that it is most appropriate to view Eq. (j45p as implementing an 
approximation to the Faddeev equation kernel in a given channel, instead of merely as a simplification 
of the propagator for the exchanged quark. From such a perspective one might argue that the Bethe- 
Salpeter and Faddeev amplitudes could influence the integrand's support within the channel under 
consideration and hence that should rightly depend on strangeness, contrary to the assumptions 
expressed in Eqs. (|58p . Such a dependence may potentially alter our expectations regarding the size of 
meson- loop corrections to the mass of baryons with different strangeness content. Notwithstanding this, 
we note that a strangeness-dependent gs is a second order effect within the rainbow-ladder framework 
because the propagator of the exchanged quark cannot itself exhibit a dependence on the environment 
in which it is subsequently embedded. Any effect should therefore be small. 

A full analysis of the impact of Eqs. (j58p can only be achieved by comparing our results with those 
obtained when Eq. (P5|) is not used. Notably, however, the pattern of results in the decuplet, reported 
here and in Ref . |6| , is consistent with the estimates of meson- loop corrections in Ref. 92| ; and hence a 
strangeness-dependent ginis not demanded. In relation to the octet some might argue, following the 
decuplet pattern or Ref. [93, that the magnitude of the meson- loop correction to an octet baryon's 
mass should fall as the number of its s-quark constituents increases. We do not find this. However, we 
can obtain to'J' — to™^ = 771^2 ^ by choosing gs = 1.15 gN- This small modification is well within 

bounds that may reasonably be applied to a first study and does not impact upon our discussion of 
band clustering, reducing masses of other S'-baryons in Table [S] by less-than 1%. 

It is also interesting to consider the baryons' diquark content, which is revealed by the Faddeev 
amplitude; i.e., the eigenvector produced by the Faddeev equation, and listed in Table[71 It is apparent 
that the relative strength of spin-zero cf. spin-one diquarks is fixed within a spectral level, almost 
independent of strangeness: the spin-zero-diquark content is 78% in ground-states; 0% in the positive- 
parity radial excitations; 52 ± 10% in the ground-state baryon parity-partners; and 46 ± 10% in radial 
excitations of the latter. This outcome is not surprising but, instead, reassuring because it indicates 
that baryon structure is largely flavour-blind within our framework, something one should also expect 
in QCD. 

This discussion highlights a striking feature of Table [T] namely, in possessing zero probability for 
J = diquark content, the Faddeev amplitudes corresponding to the first radial excitation of each 
ground-state are vastly different from the amplitudes for every other state. This suggests a fascinating 
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Table 7 Contact interaction Faddeev amplitudes for each of the octet baryons and their low-lying excitations. 
The superscript in the expression s* or a* is a diquark e num eration label associated with Eq. H3ip . except for 
[2, 3] and [6, 8], which are the 7 = combinations in Eq. H49|l . 





S' SP,3] „4 ^4 ^5 „5 ^6 ^6 „|6.81 ^|6.81 ^9 ^9 




iP = +,n = 0) N 
A 
E 


0.88 -0.38 0.27 -0.06 0.04 

0.67 -0.27 -0.45 -0.09 
0.85 -0.45 0.26 0.12 0.02 

0.91 0.14 0.08 0.39 0.00 


78% 
79% 
72% 

82% 


(P = +,n=l) N 
A 
E 


-0.02 0.52 -0.37 -0.63 0.44 

0.03 0.06 -0.78 0.63 
0.00 -0.04 0.02 0.83 -0.55 

0.00 0.01 -1.00 -0.02 0.06 


0% 
0% 
0% 
0% 


(P = -,n = 0) N 
A 
E 


0.71 -0.41 0.29 0.41 -0.29 

0.64 0.44 -0.47 0.42 
0.61 -0.47 0.23 0.55 -0.21 

0.76 -0.34 0.35 0.33 -0.28 


50% 
61% 
38% 

58% 


(P = -,n=l) N 
A 
E 


0.66 -0.41 0.29 0.45 -0.32 

0.60 0.43 -0.48 0.47 
0.57 -0.47 0.23 0.58 -0.24 

0.73 -0.34 0.37 0.33 -0.31 


44% 
55% 
33% 
54% 



new possibility for the structure of these first radial excitations of baryon dressed-quark cores, which 
was first noted in connection with the Roper resonance [93 |. 

To explain this remark, we focus first on the octet ground-states, whose Faddeev amplitudes describe 
states that are dominated by their scalar diquark components (78%). The axial- vector components are 
significantly smaller but nevertheless important. For the nucleon, this heavy weighting of the scalar 
diquark component persists in solutions obtained with more sophisticated Faddeev equation kernels 
(see, e.g.. Table 2 in Ref. [131 )• From a perspective provided by the octet states' parity partners and 
the first radial excitation of these states, in which the scalar and axial- vector diquark probabilities are 
52%-48% and 46%-54%, respectively, the scalar diquark content of the octet ground-states actually 
appears to be unnaturally large. 

One can nevertheless understand the structure of the octet ground-states. As with so much else in 
hadron physics, the composition of these flavour octet states is intimately connected with DCSB. In 
Sect. l2.2.^ we observed that in a two-color version of QCD, scalar diquarks arc Goldstone modes [85l:[86{ 
and that memories of this persist in the three-color theory. For example, in the low masses of scalar 
diquark correlations (Table HJ; and in the large values of their canonically normalised Bethe-Salpeter 
amplitudes (Table [S|) and hence strong quark-|-quark— diquark coupling within the octet ground-states. 
(A qualitatively identical effect explains the large value of the ttN coupling constant and its ana- 
logues involving other pseudoscalar-mesons and octct-baryons.) There is no commensurate enhance- 
ment mechanism associated with the axial-vector diquark correlations. Therefore the scalar diquark 
correlations dominate within octet ground-states. 

With the Faddeev equation treatment described herein, the effect on the first radial excitations is 
dramatic: orthogonality of the ground- and excited-states forces the radial excitations to be constituted 
almost entirely from axial-vector diquark correlations. It is critical to check whether this outcome 
survives with Faddeev equation kernels built from a momentum-dependent interaction. 

This brings us to a final, very significant observation; namely, the match between our computed 
level ordering and that of experiment, something which has historically been difficult for models to 
obtain (see, e.g., the discussion in Ref. |95|') a nd is not achieved in contemporary numerical simulations 
of lattice-regularised QCD (see, e.g., Ref. [9Q|). In particular, our calculations produce a parity-partner 
for each ground-state that is always more massive than its first radial excitation so that, in the nucleon 
channel, e.g., the first = state lies above the second = state. 
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A veracious expression of DCSB in the meson spectrum is critical to this success. One might ask why 
and how? It is DCSB that both ensures the dressed-quark-cores of pseudoscalar and vector mesons are 
far hghter than those of their parity partners and produces strong quark+antiquark— meson couphngs, 
which arc expressed in large values for the canonically normalised Bcthc-Salpctcr amplitudes (Tablc|31)- 
The remnants of Pauli-Giirsey symmetry described previously entail that these features are carried 
into the diquark sector: as evident in Fig. [3] and Table [S] and their comparison with Fig. [2] and Table [31 
The inflated masses but, more importantly, the suppressed values of the Bethe-Salpetcr amplitudes 
for negative-parity diquarks, in comparison with those of positive-parity diquarks, guarantee the level 
ordering we have computed: attraction in a given channel diminishes with the square of the Bethe- 
Salpeter amplitude (see App.jCj. Hence, an approach within which DCSB cannot be realised or a 
simulation whose parameters are such that the importance of DCSB is suppressed will both necessarily 
have difficulty reproducing the experimental ordering of levclsll 



5 Summary and Perspective 

We described the first DSE-based calculation of the spectrum of strange and nonstrange hadrons that 
simultaneously correlates the dressed-quark-core masses of meson and baryon ground- and excited- 
states within a single symmetry-preserving framework. The nonstrange sector was described previously 

. Therein, the physical pion mass is obtained with to„ = = m = 7 MeV; and five parameters 
are used to define the gap-, Bethe-Salpeter- and Faddeev-equations. In a comparison with relevant 
quantities, the study records a value of 13% for the overall root-mean-square-relative-error/degree-of 
freedom; and notable amongst the results is agreement between the computed masses of nonstrange 
baryon dressed-quark-cores and the bare masses employed in modern dynamical coupled-channels 
models of meson-baryon reactions. We reproduce all results in Ref. @; and in extending that work to 
hadrons with strangeness, a s-quark current-mass of 0.17 GeV produces the physical kaon mass, and 
we neither introduce new model parameters nor change the values of those determined previously. 

In explaining mesons we capitalised on recent progress in understanding the far-reaching effects 
of dynamical chiral symmetry breaking (DCSB) within the Bethe-Salpeter kernel to improve upon 
the rainbow-ladder truncation in the scalar and axial-vector channels. This enabled us to produce a 
spectrum of meson dressed-quark-core masses in which empirical values are typically overestimated 
by ^ 10%. This is a marked success since our DCSB-corrected kernel omits resonant contributions; 

1. e., effects that may phenomenologically be associated with a meson cloud. Indeed, since meson cloud 
contributions typically induce a material reduction in hadron dressed-quark-core masses, any approach 
that omits such effects, whether deliberately or inadvertently, should be viewed with caution unless 
it overestimates masses by a similar magnitude. Failing this, extreme care must be exercised before 
drawing insights or claiming understanding therefrom. 

Our description of meson dressed-quark-cores enabled a prediction of the spectrum of baryons, 
including those with strangeness, that is similarly successful. Deconstructing the masses, we arrived at 
numerous insights into baryon structure. For example, we found that the diquark content of baryons 
is largely independent of strangeness; namely, that baryon structure is flavour-blind within our frame- 
work, something one should also expect in QCD. Another noteworthy result concerns the first radial 
excitation of each ground-state; viz., they possess negligible probability for J = diquark content. 
Thus, in striking contrast to the Faddeev amplitudes of every other state, the radial excitations are 
constituted almost entirely from axial-vector diquark correlations. This possibility was first noted in 
connection with the Roper resonance [93 | . Lastly and significantly, our computed level ordering matches 
that of experiment. In particular, the parity-partner for each ground-state is always more massive than 
its first radial excitation; i.e., the first = j state always lies above the second = state. We 
showed that a veracious expression of DCSB in the meson spectrum is critical to this outcome, so that 
an approach within which DCSB cannot be realised or a simulation whose parameters are such that 
the importance of DCSB is suppressed will both necessarily have difficulty reproducing the experimen- 
tal ordering of levels. Furthermore, our results suggest that realistic simulation parameters, including 

^ A phenomenological meson mass formula is used in Ref. ^3] as the basis for the analytic inference of baryon 
mass formulae. This transmits patterns in the meson spectrum into the baryon sector. The baryon mass formulae 
inferred thereby are consistent with the empirical level ordering. This pathway bears a qualitative resemblance 
to our finding. 
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near-empirical masses for dynamical light-quarks, must be achieved in order for lattice-regularised 
QCD to confirm the role of diquark correlations in baryon structure. 

The analysis and results described herein further strengthen the claim that a symmetry-preserving 
treatment of a vector x vector contact interaction is a useful tool for the study of phenomena charac- 
terised by probe momenta less-than the dressed-quark mass. As remarked above, whilst this interaction 
produces form factors which are too hard, even they, when interpreted carefully, can be used to draw 
valuable insights, especially concerning relationships between different hadrons. This foundation is 
therefore being used in the computation of elastic and transition form factors involving baryons with 
strangeness. 

Studies employing a symmetry-preserving regularisation of the contact interaction serve as a useful 
surrogate, exploring domains which analyses using interactions that more closely resemble those of 
QCD arc as yet unable to enter. They are critical at present in attempts to use experimental data as 
a tool for charting the nature of the quark-quark interaction at long-range; i.e., for identifying distinct 
signals of the running of the couplings and masses in QCD. 
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A Euclidean Conventions 

In our Euclidean formulation: 

4 

p - q^^Piii-, (A.i) 

1=1 

{7M>7f} = 25^-- ; 7^ = 7m ; (^tJ.v = ^[7m.7i'] ; tr [757^7^ 7p 7a] = -4ep^pa ,61234 = 1. (A. 2) 
A positive energy spinor satisfies 

u{P, s) {ij- P + M) = 0^ {ij- P + M) u{P, s) , (A. 3) 
where s = ± is the spin label. It is normalised: 

m(P, s) u{P, s) = 2M , (A.4) 

and may be expressed explicitly: 

u{P, s) = VM ~ i£ [ q- ■ P ), (A.5) 

with £ = VP^ + M2, 

l)^X- = (f]- (A.6) 

For the free-particle spinor, m(P, s) = u{P, 3)^74. 

The spinor can be used to construct a positive energy projection operator: 

s = ± 

A negative energy spinor satisfies 

v{P, s) (27 • P - M) = = (27 • P - M) v{P, s) , (A. 8) 

and possesses properties and satisfies constraints obtained via obvious analogy with u{P, s). 
A charge-conjugated Bethe-Salpeter amplitude is obtained via 

r{k-P) = C''r{-k;PfC, (A.9) 

where "T" denotes a transposing of all matrix indices and C = 7274 is the charge conjugation matrix, — — C 
We note that 

Ct7^C7 = -7M, [C,75] = 0. (A.IO) 
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In describing decuplet resonances we employ a Rarita-Schwinger spinor to represent a covariant spin-3/2 
field. The positive energy spinor is defined by the following equations: 

{ij - P + M)u^(P;r) = 0, 7pU^(P;r) = 0, P^,Uf,{P;r) ^ , (A.ll) 

where r — —3/2, —1/2, 1/2, 3/2. It is normalised: 

%(P;r')Mp(P;r) = 2M, (A.12) 

and satisfies a completeness relation 

3/2 

— Yl u,{P-r)u,{P-T) = A+{P)R,,, (A.13) 



where 



2M 

r = -3/2 



1 2 - - 1 - 



with P = —1, which is very useful in simplifying the Faddeev equation for a positive energy decuplet state. 



B Bethe-Salpeter equations 

We discussed the BSEs for ground-state charged-pseudos calar a nd vector mesons in Sect. l2.2T1 and for ground- 
state scalar and axial- vector diquark correlations in Sect. l2.2i2l Here we report formulae for other channels and 
also radial excitations in each channel. In the latter context, some remarks are in order. 



B.l Radial excitations of pscudoscalar mesons and scalar diquarks 

In quantum mechanics the radial wave function for a bound-state's first radial excitation possesses a single 
zero. A similar feature is expressed in quantum field theory: namely, in a fully covariant approach a single 
zero is usually seen in the relative-momentum dependence of the leading Tchebychev mome nt o f the dominant 
Dirac structure in the bound state amplitude for a meson's first radial excitation [S^iH^iIlJ. The existence 
of radial excitations is therefore clear evidence against the possibility that the interaction between quarks 
is momentum-independent: a bound-state amplitude that is independent of the relative momentum cannot 
exhibit a single zero. One may also express this differently; namely, if the location of the zero is at fcg, then 
a momentum-independent interaction can only produce reliable results for phenomena that probe momentum 
scales fc^ < fcg. In QCD, fcg ~ 2M^ ~ (0.5 GeV)^ [23;[33;|3l. 

In the phenomenological application of a con tact interaction, however, this difficulty has been skirted by 
mea ns of an expedient employed in Refs. [98l: I99l| : i.e., one inserts a zero by hand into the kernels expressed in 
Egs. (I14|) . This means that one identifies the BSE for a radial excitation as the form of Eq. ((8| obtained with 
Eq. ([TJ and insertion into the integrand of a factor 

(B.l) 

which forces a zero into the kernel at = 1/djr, where dj is a parameter. Plainly, the presence of this zero 
reduces the coupling in the BSE and hence increases the bound-state's mass. Although this may not be as 
transparent with a more sophisticated interaction, a qualitatively equivalent mechanism is responsible for the 
elevated values of the masses of radial excitations. 

The expedient may readily be illustrated in connection with pseudoscalar mesons. First define the function 

r"(a;(Afi^, Ml a, P^)) = ^:'"(a;(M„^ , a, P^)) - 2)'"(a;(M„^ a, P^)) , (B.2) 

2 

!D'"(a;(Mi^,Mf,a,p2)) = H ds ^— ^ dr \ e^p\-TUj{Ml, M^a, P^)] , (B.3) 
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with, as usual, !Fi^(z) — —z{d/dz)f^^{z) and fi{z) — J\(z)lz. Then the BSE for the kaon's first radial 
excitation, denoted , is obtained with the kernel 



iC = 



Aee 



^Kfe 



In Eq. ((Rla)! . 



Aee Aef 
Afe Aff 



{u{MlMla,P')) 



+ 



MuMs - a(f - a)P^ - u:{MI, M^,a, P^] 



r^(u[MlMla,P'')) 



2Mr Jo 



da 



{1 - a)Mu + aM, 



fT{oj{Ml,M't,a,P^)), 



p2 



R _ 

Aef — 



I = Mr r 

Jo 



da 



[1 - a)Mu + aMs 



TT{u;{MlMla,P^)), 



-2 I ''^ 



MuMs + (f - a)Ml + aMl 



f^{oj{MlMla,P^)). 



Afe — Afe ~ Afe 



'Kff — Aff 



Aff 



(B.4a) 

(B.4b) 
(B.4c) 
(B.4d) 
(B.4e) 

(B.5) 



These subtractions are implemented so as to ensure that the leptonic decay constant of the radially-excited 
pseudoscalar meson vanishes in the chiral limit, which is a consequence of the axial- vector Ward-Takahashi 
iden tity [33 : |3^ . The kernel for the radially-excited scalar diquark is obtained through obvious analogy with 
Eq. (1231). 



B.2 Vector mesons and axial- vector diquarks 

We emphasise that the vector-meson and axial-vector-diquark Bethe-Salpeter amplitudes only assume the 
simple form in Eq. (|17|l when rainbow-ladder truncation is employed. Even with a momentum-independent 
interaction, these amplitudes possess two Dirac covariants immediately upon inclusion of next-to-leading-order 
corrections to the quark-gluon vertex; viz., 

rr (P) = 7m Si- (P) it El- (P) + i^cj^.P.F,- (P) , 7^P^ = 0. (B.6) 

Similar observations hold for a g^D{p — ~ — q) interaction [l^; |5^: 1601] . 

In studying radial exci tati ons, there are no complications in these channels: one simply implements the 
replacement C in Eq. p9[) and then works with the BSEs and canonical normalisation conditions obtained 
therefrom. 



B.3 Axial- vector mesons and vector diquarks 



Again owing to the simplicity of the interaction and truncation, the Bethe-Salpeter amplitude for an axial- 
vector meson composed from a /- and p-quark is 



ri!AP) = isjtE}tiP)- 



Inserting Eq. (|B.7P into Eq. dS} yields the following BSE: 



o = i+3^^;(-k;)^), 



2aiR 



da 



Cr{Lo{Mg, Mf, a, z)) + {MfMg + a(l - a)z)c"i [ui{M'^g , M^a, z)) 



(B.7) 



(B.8) 
(B.9) 



Given the discussion in Sect . 12. 2721 the Bethe-Salpeter amplitude for a = 1 correlation of quarks with 
flavour /, (? is given by 

r'j-''iP) = lsltEiM,AP)- (B.IO) 

In a three flavour theory the allowed correlations are {uu}, {ltd}, {dd}, {su},{sd}, {ss} and their masses are 
determined from 

^ + l^'si-mh,}^_). (B.U) 
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The canonical normalisation conditions are 
1 



-{mlt)2 

^ is ' 



F2 



(B.12) 



{/s}, 



Once again, in studying radial excitat ions there are no complications in these channels: one simply im- 
plements the replacement C — > in Eq. ([B.QP and then works with the BSEs and canonical normalisation 
conditions obtained therefrom. 



B.4 Scalar Mesons and Pscudoscalar Diquarks 
The Bethe-Salpeter amplitude for a fg scalar meson is 

r/f(P) = iD<(P). 



(B.13) 



As with axial - vecto r mesons, corrections beyond rainbow-ladder truncation do not generate new covariants. 
Inserting Eq. ([B.ISP into Eq. ^ produces the BSE 



4qir 



^:'"(a;(A^g^ Af|, a, z)) ^ Cr{oj{Ml Mf,a, z)) 
{MfMg + a(l - a)^)^l"(a;(M^^ M|, a, z)) 



The canonical normalisation condition is: 

1 



2 fi^O+ / \ 



-{mOt)2 



(B.14) 



(B.15) 



(B.16) 



It is w orth reiterating at this point that when applied in the chiral limit to the isoscalar-scalar channel, 
Eq. (fBT4|l yields [lOOl 



2Mo . 



(B.17) 



'"^(0)0++ 

Whilst this algebraic result does not persist beyond rainbow-ladder truncation, it is a useful test of symmetry 
conservation within an application of the contact interaction and also provides for a reasonable definition of a 
single dressed-quark mass-scale w hen a more sophisticated interaction is employed. 

Given the dis cussio n in Sect . 12. 2!^ the Bethe-Salpeter amplitude for a pseudoscalar diquark is readily 
inferred from Eq. ((BTSj: 

r„^J^''l(P)= Id (P). (B.18) 

In a three-flavour theory, allowed combinations are [ud], [s?i], [sd]. The associated BSE takes the simple form 

^0+ / 2 



= 1 + (-'"Mo-) 
and the canonical normalisation condition is: 

d 



El , 



3,4 



(B.19) 



(B.20) 



[/sin 



C Ground-state Faddeev equations 

C.l Z\- resonance 

As remarked when opening Sect . 13.21 owing to its inherent simplicity the A is an ideal system with which to 
illustrate steps in the derivation of Faddeev equations that involve Dirac algebra and analysis of momentum 
integrals. One begins by observing that with a momentum-indepen den t kernel, the Faddeev amplitude cannot 
depend on relative momentum and hence for the Zi-resonance Eq. (|38|l becomes 

D.p{t, P)uf{P) = /^(P) Id u^{P) . (C.l) 

N.B. Regarding Eq. (|38p in general, one might naively suppose that isospin-one tensor diquarks could play a 
material role in the Faddeev amplitude for a ground state A. However, this notion can quickly be discarded 
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because ground-states are distinguished by containing the smallest amount of quark orbital angular momentum, 
I/, and a tensor diquark is characterised hy L > 1. 
Using Eqs. (IClfl . (|45)l . Fig.[T]is realised as 



f^iP)u^iP)^Ai^ I A±MfM-,P)f^{P)u^{P), (C.2) 
with (K ^-e + p, p2 = -mi) 

3<(^; p) = 2ir,\t.}(^)ir^t..}(-P)SMz\irpl(^) , (C.3) 



where the "2" has arisen through isospin contractions 

'4 



At this point, one post-multiplies by uf{P;r) and sums over the polarisation index to obtain, using 



Eg. lTOSl) . 

,2 ^ ,4 



(2 

which, after contracting with 5^^, yields 



(P)<^(P) = 4|i^| A±M,i^l.p)At{P)R^f,{P), (C.4) 



1 = ll tro I ^ (C.5) 



-^m4[2mzi^-ii:A--P-F3M(m|„„}^^mi-h(A:-P)')]j , (C.6) 

where ^(-?^) is the canonically-normalised axial-vector diquark Bethe-Salpeter amplitude, explained in 

Eq. p7|) . Now, with the aid of a Feynman parametrisation, the right hand side becomes 

8 g?0 ^{nu}^+ f dH f\ 1 /_ 2 2 

SM^mlml^^^J (2^)4 y„ [{£ - aP^ + a{a, Mlml^^^^,ml)]^ { ^ • ^ 

+ (is: • Pf] + mA[2mA£ ■KK P + 3Af(m^„„}^_^ + {K ■ Pf)]) (C.7) 



where 

a{a, x,y, z) = {1 — a) X + ay — a{l — a) z . (C.8) 

We employ a symmetry-preserving regularisation scheme. Hence the shift £ — >■= £' + aP is permitted, 
whereafter 0(4)-invariance entails ^' • P = so that one may set 

£■ P ^ aP'^ , K ■ P = {1- a)P'^ , £ ■ K ^ a{l - a)P'^ , (C.9) 

and therewith obtain 



gio Elu}^+ r dH' f\ +(l-Q)W^)(am^ + M„) 



f'~T~^T^ f ["^{««},+ + (1 - ")^'"^][""^^ + A/^,m{„„} ,mi)) . (C.ll) 



This is an eigenvalue problem whose solution yields the mass for the dressed-quark-core of the Zi-resonance. If 
one sets gio = 1, then rriA = 1.60 GeV. 

For the purpose of illustration, here we present the Faddeev equation for the parity partner of the ground- 
state Zi-resonance: 



^ = 1^ 2""^' TTl [ -^^«]Cr(f^(«,Af^,m{ui,},_,fwi-))- (C.12) 
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C.2 Nucleon 

As we have already observed, using Eq. (|45|) the nucleoli's Faddeev amplitude simplifies to the form in Eqs. ()46|) 
which here we write as 



= s3v(P)/d, <'(P) = aivi(P)i757M + «iv2(J')75PM, i = 4,5, 



(C.13a) 



where the superscripts are diquark enumeration labels associated with Eq. (|31|l . In terms of the associated 
amplitude, the nucleon's Faddeev equation takes the form 





/ 




3Cf4i 


A.14i 


— \/2 KiiA^ 


A.142 


\ 


/ s^'iP) \ 











\/2 ^^4^141 





\/2 'KAii2 




a%iiP) 






<1 


\/2 "Ka^Ai 


%^4i 


V2 <42 


^4,2 




a%iiP) 


a%2{P) 




-^2 ^iCf^i 










^2 ^42 42 




a%2{P) 


\a%2{P)) 


V 


'KA2I 




A.4241 


■^2 3C4^42 


A.4242 


/ 


\c^%2iP)J 



(C.14) 



where the kernel expresses the isospin symmetry we have assumed; i.e., that the Bethe-Salpeter amplitudes for 
the {uu}i+ and {ud}i+ correlations are identical. Define now 



Cat 



t]s! = {1 — a) M'^ + am^ — a(l — a)m^ , 



(C.15) 



where i = 1,4: are diquark labels associated with Eq. H3ip . then the entries in the nucleon's Faddeev kernel may 
be expressed as follows: 



All 



K 



Al42 



K 



K 



A42I 



42 E 



Kee + Kef + Kpp , 



A'™ — 



/■I _ 

CnEI I daCiicrlf){amN + Mu) , 
Jo 



TUN 



A-fiV = -2cnEiFij^ I daCi(air)(l-a)(am]v + M„), 



Kpp = Cn 



2 



Al4i — -'i£;4i + J^FAi , 

E1E4 — 



"13 Jo 



da Ci(o-^)(m|(3Afu + am at) + 2Qf(l — afm%j) , 



FiEatun 



i^^F4i = ^CJV 

-"-B42 + -"F42 ) 



ml M, 



^ / da Ci(cr;^)(l - a){ml{Mu 
lu Jo 



+ SamAf) + 2(1 — a)^M„m^) : 



E42 = cjv- 



1A2 + 

EiEa , -:; 



da Ci(cr^)(Qmjv — Mu)((l — a)'^m% — ml) , 



PA2 — Civ- 



FiEa mN , t: / 4 



ml M„ 



/ da Ci{(r%){l — a){amN — — afm% ~ ml) : 

Jo 



A4il — ^AiE + -"4iF ) 



R'" — 

-114, K — 



cjv EiEa , — I 1 



3 ml Jq 



da Ci(crjv)(amjv + Mu){2mA + (1 — a) mjv) , 



JV" — 

-"-4iF — 



Cjv i^ii?4 mjv , — I 1 



3 ml Mu Jo 



da Ci(o-jv)(l — a)(2m4 + (1 — a) m,^){amN + Mu) ; 



^42B + -"■42F ) 



Cjv EiEa , , 1 



3 ml Jq 



da Ci{u^)(amN + M„)(m4 — 4(1 — a) mj^), 



J^A„ F — 



^4i4i 



A4i42 



Cjv FiEa mN , ^ , 1 



3 ml MuJo 



Cn -E4 



da Ci(crjv)(l — a)(5m4 — 2(1 — a) mN)(amN + Mu) ; 



/ da Ci(cT%){2ml(Mu — amjv) + (1 — a)^m%{Mu + 5amjv)) ; 
Jo 



■J ^AJO 



a) mN)iamN - Mu) ; 



(C.16a) 
(C.16b) 

(C.16c) 

(C.16d) 
(C.16e) 
(C.16f) 

(C.16g) 
(C.16h) 
(C.16i) 

(C.16j) 
(C.16k) 
(C.161) 

(C.16m) 
(C.16n) 
(C.I60) 

(C.16p) 

(C.16q) 

(C.16r) 



26 



3 ml Jo 



^4241 =-^^/ da Ci{a%){mi{namN+Mu)- 2(1 -aYmifi7amN + 2Mu)); (C.16s) 



3C4^42 = -^SILE^ I da Ci(cr%){m'i - (1 - afm%){. 



(C.16t) 



with El, Fi, E4, being cano nica lly normalised Bethe-Salpeter amplitudes for diquarks corresponding to enumer- 
ation la bels i = 1, 4 in Eq. (|31|) . This kernel was computed following the pr oced ure detailed for the Z\-resonance 
in App. lC.l] During this process we employed the replacements in Eq. (|C.9|) . their analogues involving the 
scalar-diquark's momentum, Jfj^^jj^^ , and K^^^-^^j^ ■ Kf^^^-^^^ — > (1 — a)^P^. In the present context, of course, 
= —m%. If one sets ~ 1, then mjv = 1.27. 

Given the structure of the kernel, it is not surprising that the eigenvectors exhibit the pattern 

"ATj = -^2a^j, j = 1,2. (C.17) 
For example, at the mass presented in Table [6l the nucleon's unit-normalised Faddeev amplitude is 



0.88 0.38 -0.27 -0.065 0.046 ' 
The axial- vector-diquark correlation provides 22% of the unit normalisation. 



(C.18) 



C.3 Kernel for the A baryon 

Here we make explicit each entry in the Faddeev equation kernel for the yl-baryon, Eq. H56|) . Define 



/ yS f,i / 2 2 N 

Ca = 4^2 jvf^ ' '^A = cr(a,Mj,m„m^) , 



(C.19) 



where aja, x, y, z) was introduced in Eq. (|C.8[) . / labels a quark flavour and i is the diquark enumeration label 
in Eq. so that mi is the mass of the associated correlation, then 



A,12 



A.161 



A.162 



!^2l 



q^A 
A.23 



A,28i 



A.61I 



n^A 
A62I 



n^A 
A.6i3 



Ca 



2MuMr Jo 
c^Ec 
M, 



da Ci{(Ja )[EiM^ - FimA^l - a)][2E2MR - FimA^l - a)][Mu + am,i] , (C.20a) 



ExMu[?,Muml + mA{ml + 2m\{l - af)a] 



F^mAil — a)(M„[mg + 2m\(l — q)^] + SmlmAa) 



Muml Jo 

Ca 
2MuMr 

Ca r 

^MlJo 

> 

2MRml 



(C.20b) 

da'Ci{cj''/)[EiMu + FimA{l - a)][ml - m\(\ - af][Mu - amx] , (C.20c) 

/ daCi{a'/)[EiMu - Fim^(l - a)][2E2MR - f2m^(l - q)][M, + amA] , (C.20d) 
Jo 



da Ci{(jY)[2E2Mr - F2mA{l - a)] [M„ -f amA\ 



c\E^ 



f 

Jo 



da Ci {cr'^'^) 



2E2MR[3Muml+mA[ml + 2ml(l - af]a] , 



-F2"i/i(l — a)[Mu{me + 2myi(l — a) ) -f Smgm^a] 



A.282 — 



c'aEg 



282 = TTT-r 2 / da Ci K^^) [2^2 A/h + F2m^(l - a)][ml - m^(l - af][Mu - amA] , 

ZMRmg Jo 

daCi{a'x^)[EiMu - FimA(l - a)][2ml + m^(l - af][Ms + amA] , 

MumlJo 



clE>, 
3M„m| 



c'aEq 



f 

Jo 

'f 

i Jo 



da Ci(cr^'^) 



EiM4ml-Am\{l-aY] 



+ Fimyi[5m6 — 2m^(l — a) ](1 — a) 



[Mu + amA] : 



da Ciial')[2E2MR - f2m^(l - a)][2ml + m\{l - aY][Mu + amA] , 



(C.20e) 

(C.20f) 
(C.20g) 
(C.20h) 

(C.20i) 
(C.20j) 
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A.62 3 



A.6181 



A.6182 



n^A 
A.62S1 



n^A 
A.6282 



i Jo 



2E2MR[ml ~ '^m\{l - af] 



+ -F2myi[5mg — 2myi(l — a) ](1 — a) 



(C.20k) 



3mf 



da (:i(cr;;''^)[M„m^[4mi - ml(l - a)^] + mi[mg + 2ml (1 - a) ](1 - a) a] ,(C.201) 



3m; 
3m| 



7' 

^'^ 422 24 4 

6 -'o 
Jo 



Mu[rn^ — 4m|m^(l — a)^ + 6m^(l ~ a)^] 



+ mA[— 9mg + 10mim^(l — a)^ + 2m^(l — a)'^]c 



^ 4 / da Ci(o-^''^)[5m6 - 7mgm^(l a)^ + 2m\{l - a)\Mu - amyi] , 
3mg Jq 



(C.20m) 

(C.20n) 
(C.2O0) 



with Ma defined in connection with Eq. and {Ei,Fi\i — 1,2}, Ee being canonically normahsed Bethe- 
Salpeter amplitudes for diquarks corres pond ing to enumeration lab els i = 1, 2, 6 in Eq. (|31|) . This kernel was 
computed following the pattern in App JC.l1 using analogues of Eq. (|C.9[) . 



C.4 U Baryon 

Given the flavour structure of the S, one can obtain its Faddeev equation directly from that of the proton by 
simply making the replacement d — >■ u. However, we assumed isospin symmetry in writing the proton's Faddeev 
equation, so that replacement is not readily achieved. We therefore report the complete structure here. 
The Faddeev amplitude for the S'^ baryon is expressed in terms of 

5^ = si (P)Id , .^f = asi [P)il^l^ + aE2 {Ph^Pf^ ,i = ^,6 , (C.21) 
where the superscripts are diquark enumeration labels, Eq. pi[l . The associated Faddeev equation is 



a%AP) 
a%iiP) 
a%2iP) 
\o^2{P)) 



A.22 ~V^A.24i A26i — V^A.242 A.262 \ 

~\/'i'K.Ai2 y/'^'KAi&i ^2 3Cfi62 

A6i2 V^A.6i4i A6161 V^A,6i42 A6162 

~\/2 ^^422 -^2 i^4'26i -^2 3C4'262 

V ^6^2 \/2 ^6^261 \/2 ^^6^42 ^6262 / 



a%^{P) 
a%r{P) 
a%2iP) 
\a%2{P)J 



(C.22) 



In order to make each entry explicit, we define 



a{a,M^l,m1,m%) 



f 



(C.23) 



where cr{a, x, y, z) was introduced in Eq. (|C.8|1 . f — u, s, i — 2, 4, 6 is the diquark enumeration label in Eq. ()31|) . 
so that rrii is the mass of the associated correlation; and then 



A.22 



K 



Ki 



K. 



24,1 F 



K22EE + K22EF + -^''22i 
-1 



K22EE = CsE2 I da Ci{crs )[Mu + ams] 
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Mr Jo 
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R Jo 



2 ml 
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24iB — Cjj- 



da Ci(crj4'')[3Msm4 + ami:(m4 + 2m%{l — q)^)] . 



''^'^ ml 2M, 



rIo 



daCiip-y, )(1 — Q:)[A'/s[m4 + 2mi'(l — q) ] + 3am4] , 



(C.24a) 
(C.24b) 

(C.24c) 

(C.24d) 
(C.24e) 
(C.24f) 

(C.24g) 
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K. 



242 E 



K. 



26iE 



K. 



26iF 



^262 



262 E 



K. 



262F 



A.4i2 



Ai2E 



K. 



^422 



K. 



422E 



K. 



42 2F 



A.4i6i 



A4261 



A.4262 



K242E + -?fM2F , 



da Ci(crJ'')[ml:(l — a)^ — ml][ams — Ms 



242F — Cjj- 



F2E4, rriE 



Kg, E + 



f 

Jo 



da Ci(cr^'')(l — a)[m%{l — a)^ — ml][amE — Ms] , 



26iB ~r -"-261^ 
E2E(i 



_ s F2-E6 TTlg 
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s E2EQ 1 ^ 

Ce 2 
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Jo 
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Jo 
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4 Jo 



4i2F — ~Ce' 



3m: 
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3mi 2Mr 
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+ amsl , 
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Jo 



jj EiEs 
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-^4i62 — 



EaEq 
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„ EjEe 
3m1m% 



f da Ci((7^'') 
io 
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Jo 
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Jo I 

+ ams{m\[—9m% + 2m%{l — af] + 2m%{l — af[Aml + m%{l ~ a)^]) 



i ^0 



mg — m|;(l — a)^][5m4 — 2m^(l — q)^][A/„ — ams] , 



(C.24h) 
(C.24i) 

(C.24j) 
(C.24k) 
(C.241) 
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(C.24n) 
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(C.24p) 
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(C.24r) 

(C.24s) 
(C.24t) 
(C.24u) 

(C.24v) 



(C.24w) 
(C.24x) 

O'f]) 

, (C.24y) 
(C.24z) 



^6, 2 — -^6, 2B + Kg 



^6i2B -fi6i2F ) 



T^E _ s E2E6 



da CiK'')[2m^ + m|(l - af][Mu 

Jo 



+ ams] : 
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tJo 
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W--"^ 3mi 2Mn. 

^622 = Kq^2E + Kq.-^2F : 
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+ amsl 



T^E _ s E2E6 
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(C.25a) 
(C.25b) 

(C.25c) 
(C.25d) 
(C.25e) 

(C.25f) 
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A.6i4i 



A,6i6i 



+ a(l — a) ms\2mg — + 2m^(l — a) ] 



Cs ^'^2^2 [ Ci(o-5'*)[m4 - m|.(l - a)^][2ml + m|(l - a)'^][ams - Ms] 



Msiml[ml + 4m|(l - af] + 2ml{l - aY[-Aml + 3mi(l - af]) 



= Cs 



EiEa 



' 3m; 

p2 /.I 



'^s n 2 



I dQCl((Tj'*) 

io 

+ ams (ml [—97716 + 8ml;(l — af] + 2m%{l — a)'^[ml + msil — ct)^]) 

/ da Ci(o-J'')[m4 — m|.(l — a)^][5ml ~ 2m%{l — a)'^][Ms — aMs] , 
Jo 

daCi{al''')[M42ml + m%(l - af) + anis{-2ml + ^m%{l - af )] , 



^6162 — '^s 



u '2E(i f 
^ 3m| J a 



da Ci(o-^' )[m6 - ms(l - a) ][ams - Mu] , 



^6261 = -f^s^-^ / rfo Ci(o-^''')[M„(m6 - 4m%{l - af) + ams{llml - 14m|.(l 
^^fi6, = c^^j— I / da Ci(o-^''')[m6 - m|(l - af][Mu - ams] ■ 



(C.25g) 
(C.25h) 



(C.25i) 
(C.25j) 
(C.25k) 
(C.251) 
(C.25m) 
(C.25n) 



with A4r defined in connection with Eq. ^ and E2, F2, E4, being canoni cally normalised Bethe-Salpeter 
amplitudes for diquarks corr espon ding to enumeration la bels i = 2, 4, 6 in Eq. H31|) . This kernel was computed 
following the pattern in App. IC.l1 using analogues of Eq. (|C.9|I . 



C.5 ^ Baryon 

As apparent from Eqs. p5|l . the baryon may be obtained from the S'^ by making the replacements u <r^ s. 
It follows that the Faddeev equation fo r the can si mply b e obta ined from that for the S'^ by making the 
replacements u ^ s and 4 -> 9 in Eqs. ((0211) . (fa22)) . ((024)) . ((OMJ- 



C.6 S* Baryon 

The Faddeev amplitude for the decuplet S* may be expressed as 



(026) 



so that the Faddeev equation can be written 



y/2%S>,\Uh{P) 



^2 3Ci 



3C 



E* 



(027) 



Defining 



■- (1 - a) M'f +ami - a{l - a)m^. , 



(028) 



^* 47r2M/ ' 

where cr{a, x, y, z) was introduced in Eq. (|08[) . / = it, s, i = 4, 6 is the diquark enumeration label in Eq. ()31|) . 
so that rai is the mass of the associated correlation; then 



E* u E^Eq 

^46 = C^. 77-2—2 



I dQCi(cr^f) 
JO 



m4(5ml + 3m^. (1 — a)^) 



+ (1 - afm%, {3ml + m%, (1 - af) 



[Mu + aniE* 



(C.29a) 



6m|m^ Jo 



m4(5mg + 3m£* (1 — a) ) 
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+ (1 — a) TTijj. (3me + m^.(l — a) ) 



S -^6 
^66 = ^E* 2 



/' 

Jo 



[Ms + ams' 



da Cii^s*) [^6 + m%,{l — af][Mu + ams*] ■ 



(C.29b) 



(C.29c) 



with £4, being cano nic aiiy normalised Bethe-Salpeter amplitudes for diquarks corre spon ding to enumeration 
labe ls i = 4,6 in Eq. (|31l) . This kernel was computed following the pattern in App. lC.ll using analogues of 
Eq. ([CH). 

The eigenvectors corresponding to the results in Table [S] are 










P = +,n = 
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P = -,71 


= 


P = = 




0.61 




0.76 




0.67 




0.66 


/I* 


0.79 




0.65 




0.74 




0.75 



(C.30) 



indicating that these states generally favour the {us}u configuration over {uu}s. This owes to the quark 
exchange character of the interaction in Fig.[T] both {uu}s and {ua}u feed into {us}u but only {us}u can feed 
{uu}s. The radial excitation of the positive-parity E* is an exception. Its origin is dynamical, connected with 
the mass ordering Mu < Ms < rrn < me, with the switch in probability occurring for 1/d.j = 0.64(2M^), at 
which point the baryon's mass is 1.85 GeV. 



C.7 S* Baryon 

As apparent from Eqs. H35p . the S*^ baryon may be obtained from the E*^ by making the replacements u s. 
It follows that the Faddeev equation for the S"" can si mply be obtained from that for the E*'^ by making the 
replacements u -o- s and 4 — >■ 9 in the equations of App. lC.6l 
The eigenvectors corresponding to the results in Table [6] are 





p 


= +,n = 





P =+,?! = 


1 


P = -,n 


= 


P = -,n = 






0.85 




0.90 




0.87 




0.87 






0.52 




0.43 




0.49 




0.50 



indicating that these states favour the {us}u configuration over {ss}u. This owes again to the quark exchange 
character of the interaction in Fig.[T] both {ss}it and {its}s feed into {«s}s but only {its}s can feed {ss}it. 



C.8 [2 Baryon 
We define 

2 

cn = j^ij^ , on := (1 - a) Aff + ami - a{l - , (C.32) 

so that the Faddeev equation for the f2 can be written 

l = 2cfi:^ r da[ml + {l-afml][amn+Ms]Ci{al). (C.33) 
rriQ Jq 

Of course, this is simply the Faddeev equation for the Zi^-resonance but for the replacement u — >■ s or, 
equivalently, 4 — >■ 9. 
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